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Abstract 

We present the complete 2-loop renormalization group equations of the supersymmetric 
standard model. We thus explicitly include the full set of -R-parity violating couplings, 
including KiLiH^. We use these equations to do a first study of (a) gauge coupling uni- 
fication, (b) bottom-tau unification, (c) the fixpoint structure of the top quark Yukawa 
coupling, and (d) two-loop bounds from perturbative unification. We find significant 

Oh, 

shifts which can be larger than the effect from the top quark Yukawa coupling. The 
value of as(Mz) can change by ±5%. The tan/3 region for bottom-tau unification and 
for the top quark IR quasi fixed point structure is significantly increased. For heavy 
scalar fermion masses 0(1 TeV) the limits on the AL ^ operators from perturbative 
unification are competitive with the indirect laboratory bounds. The two-loop correction 
to the bound on the AB ^ opertor is +15%. 



X 
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1 Introduction 

The most compelling indication for supersymmetry is the unification of the gauge coupling 
constants. This has been thoroughly investigated in the literature 0-0, mainly in the context 
of the minimal supersymmetric standard model (MSSM). In all cases it was assumed that 
R-parity 
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is conserved.^] However, R p is imposed ad hoc: there is no experimental or theoretical jus- 
tification. R p is sufficient but not necessary to guarantee the stability of the proton. Other 
discrete || or gauge symmetries are equally possible || [7| . Strict cosmological bounds based 
on GUT-scale baryogenesis have been proposed [0] but these have been shown to be strongly 
model dependent ||. The best motivation for R-parity is that it offers a good candidate 
for as yet undetected cold dark matter. However, since this has not directly been observed 
R-parity violation (ftp) should be considered on equal footing with conserved R-parity in 



super symmetric model building and especially in collider searches [10, 11 



Throughout this paper we allow for R p . Thus we consider the additional terms in the 
superpotential 

Wfe = XijkLiLjE k + X' ijk LiQjD k + \" jk UiDjD k + KiLiH 2 . (1.2) 

We have used superfield notation and the fields have the Gsm = SU(3) C x SU(2) L x U(l)y 
quantum numbers 

L: (1,2, -±), E : (1,1,1), Q:(3,2,^), £/:(3,l,|), 

D: (3,1,-|), H x : (1,2,~), H % :{1,2,\). (1.3) 

k are generation indices. Ay& is anti-symmetric in the first two indices and X"j k is anti- 
symmetric in the last two. 

Models for unification are typically constructed at very high energies, such as grand unified 
theories (GUTs). Such models can predict the absolute or relative size of parameters at the 
unification scale. In order to compare the predictions from such models with the low-energy 
data we must employ the renormalization group (RG). Recently there have been several 
studies of the RG properties of the /? p - Yukawa couplings at one- loop [112] -fll?]. These have 
been used to place first (weak) bounds on several of the higher generation operators via 
unitarity constraints [[Uj-[|l5|]. They have also been used to compare R P -GXJT predictions 



with the low-energy data [16, 17]. 

The main interest in the evolution of the gauge coupling constants is whether they unify 
at a high scale, e.g. grand-unified or string. The R p - Yukawa couplings have to date been 
neglected in this context, mainly because the effect was expected to be small. However, as we 
show the effect can be larger than the contribution from the top quark Yukawa, At, which can 
not be neglected. The main reason is that the higher generation couplings are only weakly 
bound if at all |14], |13|, |18| . Also, the bounds as presented are usually scaled with the mass 
dependent factor (rh/100 GeV) and for SUSY masses of order ITeV the bound is typically 
weaker than the bound from perturbative unification. The strictest bounds for a mass of 
0.1(1) TeV for the couplings that we will be considering are 

A323 < 0.09 (0.9), A 333 < 0.45 (4.5)*, A 3 ' 23 < 1.25. (1.4) 



1 B: Baryon number, L: Lepton number, S: Spin. 
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At 1 TeV the bound on A323 is almost identical to the perturbative limit obtained below 



in Section 3. The bound on A 333 [|18| at 1 TeV, indicated by an asterisk is obtained by scaling 
and as such is meaningless since perturbation theory breaks down for such large values. The 
appropriate bound is thus the perturbative limit. The bound on A323 is the bound from 
perturbative unification [FJ, D|. We shall thus explore all three couplings to the pertrubative 
limit. 

We propose to investigate the effect of R p on the unification of the gauge coupling con- 
stants. The Yukawa couplings only enter at the 2-loop level. We thus extend previous work 
and present the full two- loop renormalization group equations for ft p in Section 2. We then 
focus on the applications of these equations in Section 3. We investigate the unification of 
the gauge couplings, bottom-tau unification, and the fixed-point structure of \ t . We shall 
perform three case studies, one dominant operator from LLE, LQD, and UDD respectively, 
to illustrate the possible ft p effects. Given the laboratory bounds ( |1.4| ), the third generation 
couplings can be the largest and we shall focus on them for maximum effect. However, it is 
not necessarily the theoretical expectation that the third generation couplings dominate. If 
the Standard Model Yukawas and the ^? p - Yukawas have a common physical origin, a symme- 
try, then we expect this symmetry to distinguish between a Higgs superfield and the lepton 
doublet superfields. If the hierarchy in the Yukawa couplings is determined by the structure 
of the broken symmetry we would expect the resulting hierarchy to have a different flavour 
structure for the Higgs couplings then for the purely matter-field couplings || 

When studying the RG evolution of the gauge couplings the central interest is unification. 
Does it make sense to discuss R p in the context of unification? In order to avoid rapid 
proton decay we must require a symmetry which treats quarks and leptons differently. This 
is counter to any expectation from GUTs where they are in common multiplets. All the 
same, several supersymmetric grand unified models have been constructed with a low-energy 



R-parity violating superpotential [19, [2(J, S H\- These models transfer the extreme mass 



splitting in the Higgs sector to an asymmetry in the quark-lepton multiplet via LH 2 mixing. 
In order to keep this mixing small it is generated through a broken symmetry. They have no 
further fine-tuning and the proton decay rate is consistent with experiment. It is thus highly 
relevant to consider the effects of ft p on gauge coupling unification. 

In string theories there is no preference for grand unification, and unification may very 



well be obtained with a non-simple group such as Gsm- As discussed in detail in Ref.pl 
there is then also no preference for R p over for example baryon parity B p = (— ) 3B + 25 . The 
authors obtain several models with either R p or ftp. The main question in string unification 
is whether the correct value of sin 2 9w can be obtained. Again, in this context the RGEs for 
the ft p - Yukawa couplings must be considered. 

In GUTs the prediction m^Mu) = m T (Mjj) has been very successful [^, In the 
MSSM if one requires the Yukawa couplings to unify this greatly reduces the allowed region 
of the (supersymmetric) parameters. In particular one obtains a strict relation between the 
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running top mass m t (m t ) and the ratio of the vacuum expectation values tan (3 [|T], [J. Given 
the observed top quark mass |24] this results in a prediction for tan/?. How general is this 
prediction? The third generation |? p -couplings enter the evolution of m t , mb, and m T at one 
loop and can thus have a large effect. Thus if we allow for we expect the strict predictions 
of the MSSM to be modified. In Section 3.2 we shall analyze this affect and determine a new 
tan p solution for bottom-tau unification. 

With the recent discovery of the top quark we have determined all the fermion masses in 
the Standard Model. It seems that the top quark mass is special. There has recently been 
much work to predict the fermion masses at the weak scale from a simple symmetry structure 



at the unification scale p5| , B3| . It is possible that the fermion mass structure is determined 
by a broken symmetry where only the top-quark Yukawa coupling is allowed by the 
symmetry at tree-level. It's value is put in by hand, presumably of order one. The other 
couplings are then determined dynamically through the symmetry breaking model. Given 
such a model, we would then still require a prediction for the top-quark Yukawa coupling. 
An intriguing possibility is that this Yukawa coupling is given by an infra-red (quasi) fix 
point | 27[| . The low-energy value then depends only very weakly on the high-energy initial 
value; the exact opposite of a fine-tuning problem. In supersymmetric GUTs with bottom-tau 
unification one typically requires large values of At ~ 1 close to the IR quasi fixed-point. This 
has been studied in detail in Refs.[|23], [I], ||, [36]. In the MSSM this is only attained for a 
small range of parameters, in particular tan (3. We investigate the effect of the ^-couplings 
on this scenario in Section 3.4. Similar to the case of bottom-tau unification, in Section 3.3 
we find a new tan (3 solution with fixed-point structure. 



2 Renormalization Group Equations 

We apply the work of Martin and Vaughn (MV) |29| to the superpotential (|2.2|) . We shall 
closely follow their notation. At this point we retain the coupling KiLiH 2 to be most general 
and defer a discussion to Section 2.4. We denote the Gsm gauge couplings by 

#3, 92, 9i- (2.1) 



The chiral superfields are given in (|1.3| ). In Appendix A we have collected several useful 
group theoretical formulas pertaining to Gsm and the above field content. Here we mention 
that for U(l)y we use the normalization as in grand unified theories and thus use g\ — > gy. 
See the appendix for more details. We define our notation for the Yukawa couplings via the 
superpotential including all ljl v terms. 

W = (Y /,•),,/., 7/ ; Ej + (Y D )ijQiHiDj + (Y r ), ; g,//,r ; 

+ (A E k)ijLiLjE k + (A D k)ijLiQjD k + {X Ul )j k U i DjD k 

+liHiH 2 + ^UH 2 . (2.2) 
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We have introduced the twelve 3x3 matrices 



Ye, Yd, Yjj, A E k, A D k, Am, 



for all the Yukawa couplings. This implies the following conventions 



(2.3) 





^LiE}- Lj 




= (A £fc 


hj ' 




(2.4) 




Lj Ef- Li 


^E^Lj Li 


= (A £fc 


'i* ~~ 


(Aj5*)y, 


(2.5) 


Yu.DjDk 


= Y Dj fj tDk 


= Y DjDk u t 


= (Ac/, 


)jk > 




(2.6) 






— Y DkDj u z 


= (Ac/, 


)kj = ' 


"(Ac/Oifc- 


(2.7) 



We now in turn study the dimensionless couplings and then briefly also discuss the mass terms 
li, Ki. We do not here consider the soft-breaking terms. 



2.1 Gauge Couplings 

The renormalization group equations for the gauge couplings are given by 

d - 3 



_,, : _gg_£(l) + 9 1 
dt Ua 167T 2 a fl67T 2 ) 2 



i2 B abd- E (c:Tr(YlY x ) + A^TriAlA^)) 

b=l x=u,d,e \ i=l / 



:2.8i 



The coefficients B a , B ah , and C* have been given previously |30| and for completeness we 
present them in the appendix. The ^-effects on the running of the gauge couplings are new, 
we obtain 

/ 12/5 14/5 9/5 \ 
A u /> e = 6 1. (2.9) 
V 4 4 ) 

This completes the equations for the running of the gauge coupling constants at two-loop. 



2.2 Yukawa Couplings 

In general the renormalization group equations for the Yukawa couplings are given by [p£ 



d_ 

dt 



yijk _ yijp 



167T 



7 (D* + 

2 'P 1 



(16vr 



2\2^P 



(a)* 



+ (k i) + (k <-> j) 



and the one- and two-loop anomalous dimensions are 
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(i)i 



1 
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(2)j 



— 2^ mn ^ pq Y P q r Y J + y^p q y jpg g [2C a (p) — C a (i)] 



g 2 a C a (i)S a (R) + 2 2 0,(i)C 6 (i) - 3^C (i)C(G a ) 



(2.10) 



(2.11) 



(2.12) 



See the appendix for the definition of the group theoretical quantities C a (f), S a (R), and 
C(G a ). We now first give the explicit version of Eq. ( |2.10|) for the matrices ( |2.3|) in terms of 

f 1) f (2)/"- 

the anomalous dimensions, and then we present the explicit forms for 7V 3 , and 7V 3 . 

2.2.1 RG-Equations 

The RGEs Yukawa couplings are given by 

^( Y s)ifc = (Y^) ip r|* + (Y* E ) ik Tx\ + (A Ek ) ip T^ + (Y* E ) pk T L ' p , (2.13) 



d 



^(Y^)i fc = (Y^) ip r^ + (Y^) ifc r^ + (Y^) pfc r^, (2.15) 

— (A* Ek )ij = (A^p)yr^ + (A Ek ) ip T L 3 p + (y E )ikP L H x + (-h* Ek ) pj T Lp + ( Y £)jfcr| i 1 ,(2.16) 

— (A* D k)ij = (A* DP )ijT E h p + (A* Dk ) ip TQ 3 p + (A^fc)pjT^ + (Y^^r^, (2.17) 

-jj.(A-h)jk = (Av^jp^Dp + (Af/Opfe^ + (Ay P )j fe r^v (2.18) 

At two-loop the anomalous dimensions are given by 



/l ~ 16vr 27/l (16vr 2 ^ 7/l ' 



As we discuss below in section 2.4 in many cases of interest the anomalous dimension T^ 1 
vanishes. 

2.2.2 Anomalous Dimensions 

The one-loop anomalous dimensions are given by 



(i)Lj 
Tl- 


= (Y E YU.. + (A E ,A*,,) (j 4 


3(A Eq A D q)ij 


-Si(-g Y + -g 2 2 ), 


(2.20) 




= 2(Y E Y E ).+Tr(A E ,A E3 


)-^'(^y), 




(2.21) 




= ( Y ^)« + ( Y ^) y + 


{A m ADi)ji — 


c* 7 / 2 ^2 ^ 2 \ 

Oi(^g Y + 2^2 + 3^3)) 


(2.22) 




= 2 (Y^Yo).. + 2Tr{A D ^ DJ ) + 2{A m A\ Jq ) ij ~ Hi^r + |?f)), 


(2.23) 




= 2(Y^Y [/ ).. + Tr(A t/i A t af ; 


+ 




(2.24) 




= Tr (3Y D Y^ D + Y E Yi)-( 


3 9 3 2\ 




(2.25) 




= 3Tr(Y u Yl)-(^ Y + 3 - 


gll 




(2.26) 




= 3(A D qY* D ) iq + (A E qY* E ) iq . 






(2.27) 



For the two-loop anomalous dimensions we write 

lf" = hf !> ) „ • (2-28) 

\ " / vukawa \ 1% / g— y \ J1 'gauge 



These correspond respectively to the three terms of (|2.12|) . These are given explicitly below. 
The pure gauge two-loop anomalous dimensions are given by 



/ (2)LA 

W^i J 

\ 1 / gauge 




j 15 4 207 4 9 2 2\ 

^ 4 52 100^ y k/ 2 ^' 


{ c\ c\r\\ 

(2.29) 


UEi J 

\ 1 / gauge 




« 23 4 4 
** -25 ^ 


(2.30 J 


(7§*) 

V ^ l /gauge 




ci / ^4 1^4 1 9 9 _4 o 2 2 22 22\ 

g ^3 + 4 ^2 + goo ^y + &9 3 92 + 45 9 3 9 Y + ^yy), 


(2.31) 


ceo 

\ 1 / gauge 




5H--ai + — o* + — <7 2 Q I) 
A 9 9i 225 9y 45 939y) ' 


(2.32) 


{■$"') 

\ 1 / gauge 




■ 8 4 856 4 128 22 v 
il 9^ 3 + 225^ y+ 45 539yj ' 


(2 33) 


V 1 / gauge 




(tS" 2 ) = hS )Ls ) , 

V 2 / gauge V 1 / gauge 


(2.34) 


(■>?!"•) 

\ 1 / gauge 




o, 


(2.35) 


The mixed gauge- Yukawa 


two-loop anomalous dimensions are given by 




V 1 /fl-l/ 




(16s 2 - -c? 2 ) (A D ,At fl ) + ~ 5 2 (Y E Y^ + A^A^)ii, 


(2.36) 






M - ^y)(YiY^)« + (3g 2 2 - lg Y )Tr(A E , AL), 


(2.37) 






-oEfTYnYt + 2YrrYt"l + fAt,Ar> ^ 1 


(2.38) 






(^-9 2 + ^y) (At/sAjA. 
3 15 v J v 

+ (%2 + 2 -g Y )[(Y D Y ] D ) . . + Tr(A. D iA) Dj )], 


(2.39) 






(6^ - |&) (Yf/Yt,) .. + (^ 3 2 - l^jTrCA^A^), 


(2.40) 






(16^ 2 - \g Y )Tr{Y D Y^ D ) + ?<&Tr(Y s Yi), 


(2.41J 






(16s 2 + ^-gDTr {Y v Yl) , 


(2.42) 






(16^ 2 - ^ 2 ) (A D9 Y^ g + ^ 2 (A E9 Y%) iq . 


(2.43) 


The pure Yukawa two-loop anomalous dimensions are given by 




-(72*) , = 

\ 1 / yukawa 


2 (Y 


■ E YlY E Yl) . . + (Y E Yi) . . Tr (Y b Y], + 3Y D Yy 




+ 








+ 




) im (y£a b ,A^) . + (A £n At) . . Tr (At A^) 
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+ 2 (A £ nA£ r ) i3 . (YtY B ) nr + (A E mY* E Y E A* Em ) i:j + (A^A^Y^ (Yj,) m . 

+ 3 (A E mA* Dq Y D ) iq (Y E ) m . + 3 (A^A^A^A 

+ 6 (A D nAtr) . . ((YtY D ) nr + (A^A^) nr + Tr (A^A D ,)) 

+ 3 ( A^jm Y^Y^A^m ) + 3 ( Apm Y^jY^A^m ) + (A^m A^qA^q Ag m ) 

\ / ij \ / ij 

+ 3 (ADmAS,,A Dfl A{, m ) , (2.44) 

= 2(Yt Y£ ) ji Tr(YtY E + 3Y i3 Yt ) ) + 2(YtY E YtY B ) jj 

+ 2 (Yj.A B ,A^Y B ) + 6 (YjjA^A^.Ys) + 2 (Y^A^A^Y £ ) (2.45) 



+ 6 (YtA^A^Y^) .. + 2Tr (A^Y E Y^A^) + 2 (Y^A^A^Yb)^ 

+ 6 (Y^A £i A^Y D ) , g + 2Tr (A E i A* Eq A E g A Ej ) + 6Tr (A^AcA^A^) , 

-(^U wa = 2(Y D YtY z? Yt) y . + (Y D YL)^Tr(YtY £ + 3YLY D ) 

+ 2 (YuYhYuYh) . 3 + 3 (Y„ Y*,) y Tr (y|, Y„) 

+ 2 ^YijA^-gAj/gY^ + 3 ^A^mAogA^qA^mj^ 

+ 2 (Y D ) m (Yj,) rj Tr (A^A^) + (A^Y E Y^A^) .. 

+ ( A Dm A^ 9 A Eq A D m J + (Yo) jm ( 3Y|)A^,AJ, m + Y^A^A^,™ J 

+ (Y u ) iB (Yj,) r .Tr(A^A^) 

+ 2 (A^A^n) ((Y{,Y D + A^A^) w + Tr (A^A D ,)) 

+ (Yy m .(3A^A^Y D + A5 m A^Y E ). 9 , (2.46) 

"(7S Dj ) . = 2(Y5,Y D Y{,Y D ).. + 2(Y5,Y t/ Y| / Y D ).. 

V ' / yukawa \ /ji V ' Jt 

+ 2 (Y{,Y D ) Tr (Y^Y S + 3Y^Y D ) + 2 (y dA^A^Y^) 

+ 2 (y^A^A^Yd)^ + 6 (Y^A^Y^. 



+ 2Tr(A^A^(Y D Y 

d + Yf/Yjj] + A D iA Dq Ai)qA Dj 
+ Tr (fSAflq A^ g A D i A^ 3 + 2A£>iA^ Ag 9 A^ g + 2Y B Y^A^A^) 
+ 2 (YjjA^A^Ys) . +6(Yi,A^A^Y D ). 

\ / jq \ J jq 

+ 4(A Um ) in (A\j m ). r Tr(A<< Dn A D r 

+ 4 (Aljm (Y^ D Y D + A*jj q AjJq^ A* v 



+ 2(A U nA* ur ) ij (2(Yl J Y u ) nr + Tr(Al n A U r)), (2.47) 
-b% Uj ) yukawa = 2 (YtY^Yc) „ + 2 (Y^Y^ + 6 (Y^) „ Tr ( Y „Y&) 
+ 2(Y + 4Tr (A c^jAf/iYjjYi)^ 
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+ 4Tr (A^A^A^A^) + 4 (a^A^IY (a^A d ,) , (2.48) 

+ Tr (3\ Dq \^ Dq Y E Y j E - A Eq \ j Eq Y E Y^ E - QY^ D Y D A\j q A uq 

+ ?,Y D Y^ D X ] Dq K Dq ) + (Y E Y E ) rn Tr (\ En A E r) 

+ 6(Y\,Y D ) rn Tr(A^ n A D r), (2.49) 

~hf R ) , = TrfaYuYlYuYl + ZYuYlYnY^+SYuYlAl^) 

\ 11 / yukawa \ / 

+ 3(Y^Y c/ ) TO Tr(A^A l/ ,), (2.50) 
- (7i? Hl ) . = 3 (YaYtA^YJ,) . + (Y B YkA| ? Y*) . 

\ 1 /yukawa \ / iq \ / iq 

+ 3 (YtY^YtAin)^ + (A £ .Y*). Tr (A^A^) 

+ (Afim Agq AfisY^)^ + 3 (Y E Ad<1 A^gA^m) . 

+ 9 (Yj,Y D Yl,AS„) n . + 6 (Y f D A E „) r . Tr (A^ A D r) 
+ 6 (A^A^Y^A^)^ + 3 (Y+jYcY^A^)^ 

+ 3 (A^A^A^Y^) . (2.51) 

This completes the renormalization group equations for the Yukawa couplings at two-loop. 
Before we discuss applications we briefly consider the renormalization of the bilinear terms. 

2.3 Bi-Linear Terms 

Following the general equations given in MV the renormalization group equations for the 
bilinear terms now including all R-parity violating effects are given by 

4* = /4r£ + rt}Wrt, (2.52) 



dt 
d 

—rv — rv i ff2 t nr 1 Lp T ^i- Hl 



«< = kTZI + k^+^J. (2.53) 



The anomalous dimensions at two-loop are given in the previous subsections. As already 
noted in MV the bi-linear terms do not appear in the equations for the Yukawa couplings. 

2.4 Discussion 



Equation ( |2.53j ) implies that for «j = at tree-level for all i, a non-zero Ki is generated via 



the //-term. However, as is well known, if the coefficient of the corresponding soft breaking 
term equals that of the superpotential term, then the terms KiLiHi in the superpotential 



can be rotated away through a redefinition of the Li and Hi [|19| . If we are considering 
the one-loop or two-loop renormalized Lagrangian then we must make this rotation after 
renormalization. Again, there will be no term KiLiHi in the Lagrangian, the rotation matrix 
will of course differ. Thus there is no mixing between and Hi and 1 = is guaranteed by 
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the relevant counterterms at two-loop. This then also applies to the Eqs ( p . 1 31 , 12 . 1 4| , [2 . 1 6| , |2 . 1 7\ ) 
in Section 2.2.1 and Eg. ( |2.53|) in Section 2.3. If the Lagrangian has additional symmetries 
which distingiush between Li and Hi and allows LiB.2 in W^ p then these terms must be 
retained. 

The two-loop renormalization group equations for the Yukawa couplings also respect sev- 
eral symmetries. If at some scale for example A" jfc = for all i, j, k then baryon parity, B p , 
is conserved at this scale. There are no ^-couplings in the theory and thus in perturbation 
theory no ^-couplings are generated, i.e. the RGEs preserve A^ fc = 0. Analogously lepton 
parity, once imposed, is also preserved by the RGEs. If at some scale Xijk = X"j k = for 
all i,j,k and only one lepton flavour is violated, e.g. A 3j7c ^ then this is also true for all 
scales, provided rj^ = 0. So the lepton number violating couplings LQD and LLE decouple 
in this limit. If the neutrino masses are non-zero then this is no longer true. The electron 
mass matrix Y^ then contains off-diagonal entries which generate off diagonal T^ } T^ via 
the RGEs. But the effects will be very small and can thus be neglected in most circumstances. 
However, if we assume only A' m ^ at some scale then through the quark CKM-mixing the 
other terms X' U j will be generated. 

Our results agree with MV for the MSSM Yukawa couplings. We also agree with the 



one-loop results 13, 12, 14 



3 Unification 

We now apply our two-loop RG-equations to the questions of unification. We shall assume 
as a first approximation that the ^-couplings have a similar hierarchy to the SM Yukawa 
couplings and thus only consider one coupling at a time. The third generation couplings have 
the weakest bounds ( |1.4| ) and can thus lead to the largest effects. We shall consider the three 

cases 

LLE : A 323 , LQD : A 333 , UDD : A 3 ' 23 . (3.1) 

We assume that in each case the respective operator decouples from the other |? p -operators 
whose couplings we set to zero. Is this approximation consistent, i.e. is it stable under the 
RGEs? First consider the LLE operator. The coupling A323 only violates L M and could thus 
also generate non-zero A 121- However, A121 decouples from A323 if it is the only non- vanishing 
operator at some scale. The other LLE operators are protected by global L e and L T . It is 
thus consistent to only consider a non-zero A323. 

The coupling A 333 will generate X' 3i j 7^ in higher order. However, these terms are propor- 
tional to third generation off-diagonal CKM-matrix elements which are very small and can 
be safely ignored. Thus for A 323 the decoupling assumption is also good. For A 323 the higher 
order mixing with the couplings X'- 2k is also very small. The mixing with A 313 is Cabbibo 
suppressed and involves a first generation Higgs Yukawa coupling and can be neglected. Thus 
in all three cases the decoupling is a good assumption. In line with this argument we assume 
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the following form for the Higgs- Yukawa matrices 

Ye = diag(0, 0, A r ), Y D = diag{0,0, X b ), Y v = diag(0, 0, A t ). (3.2) 
In order to determine the scale of unification we numerically solve the renormalization 



group equations. Our analysis is analogous to that of Refs |3T[] to which we refer for de- 
tails. We also do not consider GUT threshold corrections. Our analysis differs in that we 
restrict ourselves to three generations but instead add in turn one of the three ^-Yukawa 
couplings (|3.1|) . Thus we run the full set of equations including the two- loop correction of 
the Yukawa couplings to the running of the gauge couplings (|2.8j). We use the experimental 
values OL^iMz) = 127.9 and sin 2 W (M Z ) = 0.2324 0. We consider only a single super- 
symmetry mass scale Msusy — Mz- F° r more realistic spectra Msusy can be considered 
an effective mass scale which enters in the RGEs 0, [|. The supersymmetric mass spectrum 
can be highly non-uniform with the masses typically larger than Msusy- We then iteratively 
determine the value of a^(Mz) and thus also the unification mass scale Mu and the coupling 
at the unification scale a\j We have chosen a running top mass of mt(m t ) = 165 GeV 
which corresponds to a pole mass m t = 175 GeV in agreement with the discovery at the 
Tevatron We fix the running bottom quark mass at mf,(m&) = 4.25 GeV |3^| . 

We solve the RGEs for different values of the ^-coupling at the weak scale, starting from 
zero. The maximal value we consider is where the running coupling reaches the perturbative 
limit at the unification scale. So we require 

A| (Mu) 

KV ' < 1 (3.3) 
An 

or \jt p (Mu) < 3.5. 

One of the encouraging aspects of supersymmetric grand unified theories is the possibility 
of bottom-tau unification at Mu. Similar to gauge coupling unification this is not possible 
in standard GUTs. There has been a large interest in the literature |22|, [l], H [3i| in the 
restrictions on the unification scenario from bottom-tau unification. Requiring bottom-tau 
unification leads to a strict relation between the running top quark mass and tan j3. For the 
experimental value of m t [21"|, tan /3 is predicted to be very close to 1.5 or around 60 We 
are interested in how the effects of can relax this strict relation and allow a larger range of 
tan p. As a model scenario we consider tan (3 = 5 which is well away from the solutions in the 
MSSM. We then investigate the possibilities for a top fix-point solution as well as bottom-tau 
unification including the ^-effects. We do not consider the other GUT mass ratios. 



3.1 Gauge Coupling Unification 

Before determining the effects of X^ p ^ we would also like to consider the effect from the 
non-vanishing top quark Yukawa coupling. Thus we first determine the unification values for 
Ai = and = 0. We obtain 

a 3 (M z ) = 0.128, Mu = 2.98 10 16 GeV, a v = 0.043. (3.4) 
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In order to discuss the effects of the non-zero Yukawa couplings we consider the unification 
parameters as functions of the weak scale values of At and X^ p 



a s (M z , At, Xftp), M u (X t ,X ftp ), au{X t ,X^ p ). (3.5) 

and define the ratios 



Ra 3 {Xt, Xj 



a 3 (M z , At, X np 
a 3 (M z ,0,0) 



R Mu {X u X K ) - , (3.6) 



-Ray (At, Xft p/ 



&u(Xt, AgJ 
ac/(0,0) 



We first consider X$ p = and fix X t (M z ) = 0.44 so that m t (m t ) = 165 GeV (tan/5 = 5). 
This shifts the ratios (|3.6j ) as is shown in Figs la, 2a, and 3a, Ri(X t ,0) ^ 1. It can be read 
off for (A £ 3) 23 (M Z ) = 0, {A D3 ) 33 (M Z ) = 0, or (A C73 )2 3 (M Z ) = respectively. The shift is of 
course identical in the three plots, note however the different scales for Ri. Due to the large 
top Yukawa a s (M z ) and ay are lowered only by about 1%. The unification scale is lowered 
by about 1.5%. 

Next we turn on the ^-couplings. In Figures lb, 2b, and 3b we can read off the value 
of the ^-coupling at the unification scale as a function of the coupling at the weak scale. 
{A-e 3 )23{Mjj) reaches its perturbative limit (|3.3|) for a weak scale value of 



(A E s) 23 (M z ) = 0.9. (3.7) 

It is worth pointing out that this is the same as the laboratory bound for slepton masses at 
1 TeV\ Thus although the laboratory bounds on the LLE operators are generally considered 
to be very strict; for heavy supersymmetric masses they are no stricter than the perturbative 
limit. At this point (A£3) 2 3(M[/) has run off the plot but it should be clear how it extrapolates. 
All other quantities in Figs. la,b are plotted up to (A E 3) 23 (M^) = 0.9. We have chosen the 
scaling of the plot so as to highlight the effects on the other quantities. The perturbative 
limits for the other couplings are given by 

(A D3 ) 33 (M Z ) = 1.08 (3.8) 
(A U3 ) 23 (M Z ) = 1.14 (3.9) 



The latter value is about 15% higher than the 1-loop value previously obtained [|13|, [14| . This 
is the same order as the two-loop Yukawa corrections obtained in p]. 

In Figures la, 2a, 3a we show how the ratios ( |3.6|) change as we turn on the ^-couplings. 
For (A E 3) 23 (M Z ), a s and au are practically unchanged except very close to the perturbative 
limit. However, M v is shifted upwards by up to 10%. For (A E 3) 23 (M Z ) pa X t (M z ) the 
decrease due to the top quark is off-set by A# . For (A D 3) 33 (M^) the maximum combined 
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shift in a s (Mz) is a decrease of 4% giving a value of a s (Mz) = 0.123, in better agreement with 
the data [j35] . otu is decreased slightly. However, M\j is decreased by up to 20%. This effect is 
significantly beyond the effect due to the top quark Yukawa coupling. For {K u -i)2z{Mz) 1 aiu 
remains practically unchanged. a s now has an overall increase of about 5% at the perturbative 
limit corresponding to a value of a s [Mz) — 0.134 in disagreement with the experimental value. 
Mjj is raised by up to 20%. 

Thus we find ajj essentially unchanged by ^-effects. M\j can change either way by up 
to 20%. If we compare this with other effects considered in Ref. || we find it of the same 
order as the uncertainty due to the top quark Yukawa coupling or the effects of possible 
non-renormalizable operators at beyond the GUT scale. The effect is much smaller than that 
due to GUT-scale threshold corrections or weak-scale supersymmetric threshold corrections. 
It is thus much too small an effect to accomodate string unification. The strong coupling can 
also change either way by up to 5%. A decrease is favoured by the data and is welcome in 
supersymmetric unification. The effect of the .^-couplings on a s [Mz) is of the same order as 
the effects due to the top-quark Yukawa coupling, GUT-scale threshold effects and high-scale 
non-renormalizable operators 0. 



3.2 b-r Unification 

In order to study the unification of the bottom and r Yukawa couplings A&, A r at Mjj we 
define the ratio 



r Wt (Mu) = ; b ) n ; , (3.10) 



\ b (Mu,\ np ) 
For = 0, tan (3 = 5 we have 

R b/T (Mu) = 0.74. (3.11) 

Thus including the top-quark effects but before turning on the .^-coupling we are well away 
from the bottom-tau unification solution R b / T (Mu) = 1. Recall that the uncertainties due to 
the bottom quark mass are small for small tan (3. Now we consider the corrections due to the 
-^-couplings. The one-loop RGE for R b / T {t) is given by 



dRb/r(t) _ R ( s 



+ ~ 3 ^r - 3 ^323 + 2 ^323 + ^9y ~ 



(3.12) 



At one-loop the evolution of R b / T is unaffected by the operator L3Q3-D3. The slight rise which 
we observe in Fig 2b is a two- loop effect which is very small. The leading dependence of R b / T 
on A323 has a negative sign and as we see in the two- loop result shown in Fig. lb R b / T drops 
significantly. At the perturbative limit it has dropped by a factor of three and R p becomes 
a dominant effect on the evolution of R b / T . This is important for the range of tan/5 which 



leads to bottom-tau unification. In the MSSM R b i T is too large for tan/5 ~ 1.5 or ~ 60 fl3~4 



b/r 

Including a non-zero operator A323 strongly reduces R b / T and thus can lead to bottom-tau 
unification outside the previous regime. 
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For A323 7^ there is an additional positive contribution in the evolution of Rb/ T (t). The 
full two-loop result shows a clear rise in Ri/ T (t) as a function of A323 in Fig. 3b. The maximum 
increase at the perturbative limit is by 60%. For X 323 (Mz) = 1.07 bottom-tau unification is 
restored! This is quite remarkable. Even though /^-couplings are usually expected to lead 
to only small effects they can have a significant impact on our understanding of Yukawa- 
unification. Recall, that grand unification is possible in ^-theories Jl9| . 



From the structure of the couplings it should be clear that for example for X' 233 we get 
an increase in Rj,/ t (Mjj) as well since at one-loop it does not couple to A T . Similarly, A 322 
decreases Rb/ T - This leads to further bottom-tau unification points. 

3.3 Fixed Point Structure of the Top Yukawa 

We would next like to discuss the effects of the .^-couplings on X t {Mu). In particular we are 
interested in finding fix-point structures. These are defined by zero derivative at Mz- For the 
MSSM parameter-point we have chosen, tan/3 = 5, we obtain 

\ t {Mjj) = 0.45, forA #p = 0. (3.13) 

This is indicated in Figs lb, 2b, and 3b. The one-loop evolution equation of X t is given by 

(3.14) 



QXf + X 2 b + 2A 3 ' 2 2 3 + X'i i3 - (^g 2 Y + 3g 2 2 + ^-g 2 3 ) 



from which we can read off the dominant effects. At one-loop, the evolution of At is inde- 
pendent of the LLE operators. If we look at equation |2.15| , as well as the equations for the 



two-loop anomalous dimensions Tq^, T^, and rj| we see that for A- jfc = 0, for all i,j, k, there 
is also no dependence on LLE at two-loops. Thus X t (Mu) is independent of X^{.Mz) at 
two-loop, as we see in Fig. lb. 

For L 3 Q 3 D 3 we see that X t (Mjj) should grow with X 333 . The full two-loop effect is shown 
in Fig. 2b. At the perturbative limit At is large but clearly unequal to 1. The increase is by 
about 60%. The derivative at Mz does not vanish and we do not observe fix-point structure 
for At or A 333 . For lower values of tan/3 but well above 1.5 new fixed point solutions exist. 

For U^D 2 D 3 we have a larger positive coefficient than in the previous case and as we see 
in Fig 3b At rises more quickly as we approach the perturbative limit. For X 323 (Mz) ~ 1, X t 
has the fix-point structure, i.e. 

1 a 

QX 2 (M Z ) + 2X' 3 %(M Z ) » -47i 2 a 3 (M z ). (3.15) 
This is roughly the same point where R^/ T m 1. Thus for U3D2D3 bottom-tau unification 



also corresponds to top IR fix-point structure, as in the MSSM |27|, |3(| |TJ. This is all quite 
remarkable. The /? p -couplings can have significant affects on the entire Yukawa unification 
picture. From the one-loop RGE for X 323 we can see that we are well away from its infra-red 
fixed-point. 
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4 Conclusions 



We have argued that R p is theoretically on equal footing with conserved R p . Since it can 
be realized in grand unified theories it is relevant for unification. We then first determined 
the complete two-loop renormalization group equations for the dimensionless couplings of the 
unbroken supersymmetric Standard Model. It is only at two-loop that Yukawa couplings affect 
the running of the gauge coupling constants. We then considered three models of R p . We 
have added to the MSSM in turn the three Yukawa operators L 3 L 2 E 3) L 3 Q 3 D 3) and U 3 D 2 D 3 . 
We considered their effects on various aspects of the perturbative unification scenario. We 
have focused on qualitative effects. A detailed search for a preferred model is beyond the 
scope of this paper. We found several important effects. The unification scale is shifted by 
upto ±20%. This is comprable to some threshold effects but insufficient for string unification. 
a s (M z ) can be changed at most by ±5%. The reduction which is favoured by the data is 
obtained close to the perturbative limit of A 333 . We have obtained the two- loop limit from 
perturbative unification for all three operators. For A323 it is equivalent to the laboratory 
bound for a slepton mass of 1 TeV. For A323 the two-loop limit is 15% weaker than the one- 
loop limit previously obtained. For bottom-tau unification we have found significant affects. 
For L 3 L 2 E 3 bottom-tau unification can be obtained for values of tan/3 < 1.5. For U 3 D 2 D 3 
we found a new point of bottom-tau unification at tan (3 = 5. It is remarkable that this point 
also shows top Yukawa infra-red fixed point structure but no \ 323 fixed point structure. 
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Appendix 

We consider a group G with representation matrices t A = (t)f J . Then the quadratic Casimir 
C(R) of a representation R is defined by 

(t A t B )i = C(R)6i. (A.l) 
For 577(3) triplets q and for SU(2) L doublets L we have 

CW(3)(?) = |, C SU(2) (L) = ^. (A.2) 

For U(l)y we have 

C(f) = (A.3) 

where Y(f) is the hypercharge of the field /. The factor 3/5 is the grand unified normalization. 
For the adjoint representation of the group of dimension d(G) we have 

C(G)5 AB = f ACD f BCD , (A.4) 
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where f ABC are the structure constants. Specifically for the groups we investigate 



C{3U{3) C ) = 3, C(SU(2) L ) = 2, C{U{l) Y ) 
and C(SU(N)) = N. The Dynkin index is defined by 

Tr R (t A t B ) = S(R)5 AB . 
For the respective fundamental representations / we obtain 

SU(3), SU{2) : S(f) 
U(l)y : 



1 

2' 



where we have inserted the GUT normalization for U(1)y- 

The coefficients in the two-loop running of the gauge couplings 



(A.5) 
(A.6) 

(A.7) 
(A.8) 



BP 



B 



(2) 
ab 



c. 



u.d.e 



(f ,1,-3), 

/ 199/25 27/5 88/5 \ 

9/5 25 24 
V 11/5 9 14 J 

( 26/5 14/5 18/5 \ 

6 6 2. 

v 4 4 J 



are given by 

(A.9) 
(A.10) 

(A.H) 
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Figure 1: We have fixed the values tan/5 = 5, mb{rrib) = 4.25 GeV and m t {mt) = 165 GeV. 
For a vanishing top-quark Yukawa coupling and a vanishing .^-coupling we obtain the uni- 
fication values shown in Fig (la). In (la) we then show the ratios Rm V i Ra v and R a:j at 
the unification scale as a function of (A E 3) 23 (Mz). In Fig. (lb) we show the values of Rb/ T , 
At and (A E 3) 23 at the unification scale as a function of (A E 3) 23 (Mz). (A E 3) 23 (Mt/) reaches 
its perturbative limit for (A E 3) 23 (Mz) = 0.9. At this point it is outside the plotted region. 
The plotted region in Figs (a) and (b) is chosen so as to highlight the evolution of the other 
quantities. 

Figure 2: The same as Figure 1 just for (A^3) 33 except the perturbative limit is reached for 
(A DS ) 33 (M Z ) = 1.08. 

Figure 3: The same as Figure 1 just for (A l /3) 2 3 except the perturbative limit is reached for 
(A l/3 ) 23 (M Z ) = 1.14. 
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